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Abstract 

-»,. ' We set up a discrete space-time dynamical model of molecules with 

Qs I thermalised kinetic energy and repulsive cores, in an external poten- 

tial. The model obeys the first and second laws of thermodynamics. 
The continuum limit, obtained using a MAPLE program, gives rise to 
coupled reaction-diffusion equations for the density and temperature 
1 -^ I fields. The system obeys Onsager symmetry and exhibits the Soret 

C^ ' and Dufour effects. 
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1 Introduction 

It was known from experiments in the Nineteenth Century that a liquid at 
uniform temperature, but with a gradient in its concentration, develops a 
non- uniform temperature (the Dufour effect); this was confirmed by Wald- 
mann [^^. Thus Fourier's law |13] needs modification. A local version of 



the Thompson effect [26| was also found, in 1856, by Ludwig: if the density 
of a liquid is uniform, but the temperature is nonuniform, then there is a 
transport of matter. Thus, Fick's law jl^ needs modification. The ratio 
of the molecular migration to the temperature gradient is called the Soret 
coefficient, after C. Soret, who studied the phenomenon in 1879-81 [19|. For 



gases, the Soret effect, under the name 'thermal diffusion', was obtained the- 
oretically in 1911 by Enskog |ll| using kinetic theory, and also by Chapman 



in 1912 1^. In Enskog's work, it showed up in the Lorentzian gas but not 
in the Maxwelhan gas. It was first observed experimentally in gas mixtures 
by Chapman and Dootson Q. 

Concerning Enskog's work, Hirschfelder et al. say "Each time we lower 
the level of description it is necessary to introduce a condition which restricts 
the possible states under consideration. In this case [Enskog's method] it is 
not clear how the restriction has been imposed" [0], p. 492. In the present 
paper we derive a model of a dense fluid using the methods of statistical dy- 



namics 1 21], which is a systematic way to reduce the level of description. The 
possible states are restricted using information geometry in a well-defined 
way [^, 1^. In the present model, the essential part of the dynamics, the 
random hopping of molecules to neighbouring holes, is included, but there 
is no interparticle potential. The kinetic energy of the particles is fully ther- 
malised, in the spirit of Smoluchowski ^. We find that the model exhibits 
both the Dufour and the Soret effects, while obeying the first and second 
laws of thermodynamics. This success seems to contradict the statement of 
jy], p. 103: "No really satisfactory simple theory of this thermal diffusion 
can be given... The reason is that thermal diffusion is an interaction phe- 
nomenon. Similar remarks apply to the inverse 'diffusion thermo-effect' [the 
Dufour effect in gases]". Our theory, statistical dynamics, is not as simple 
as that using free paths referred to in this quotation, and it may be a mat- 
ter of opinion whether it is 'really satisfactory'; however we do show that 
the effects follow from the assumptions that the state is in local thermal 
equilibrium, and that the hopping rate is proportional to the kinetic energy. 
These are kinematic assumptions, and do not require the solving of a model 
with an explicit interaction between the particles. The interaction enters 
only implicitly; its effect is replaced by the hopping term and the exclusion 
principle, followed by local thermalisation. So the last part of the quotation 
is not true. 

In our model, the system is described by the particle density p{x, t) and 
the temperature field Q{x, t). The potential energy of a particle at x is V{x)^ 
and the heat capacity is unity. Thus the density of heat is pQ. There is a 
maximum possible density, denoted pm] this corresponds to a hard core of 
diameter £, where pm = (■~'^ in v dimensions. The particle current jc, and 
the heat current j^, are given in terms of /j, G by 

je = -\{QVp + p{l-p/pm)y{Q + V)) (1) 

j^ = 2 (ej, - Ap(i - p/p^)Gve) . (2) 

Here A is the microscopic hopping rate. These relations in conjunction with 



the conservation laws 

|^ + divje = 0; ^M + divj^ = -je.Vy (3) 

determine the dynamics. The heat current j^ is not conserved, because of 
the heat source —jc-^V = jc-F, where F is the force. Thus the work done 
by the field is entirely converted into heat; we call this the Smoluchowski 
point of view, though it is implicit in |^. The current 

je = h + ^ic (4) 

of the energy density p{@ + V) is conserved, as it obeys 

^/5(G + V)+ divj, = 0. (5) 

Thus the system obeys the first law of thermodynamics. 

The particle current jc carries with it, by convection, a heat flow of Qjc- 
This already suggests that the Dufour effect is to be expected. The surprise 
here is that the Dufour effect is 20 jc, double what is expected from this 
intuitive argument. We call the difference the 'anomalous convection'. The 
Soret effect comes from the term — A/?(l — p/pm)V0 in jc- This remains 
non-zero as pm ^>- oo, so the effect does not depend on the presence of a 
hard core. We see that in this limit, the Soret coefficient jc/(pVB) is A. 
The fact that in eq. (g) the temperature is added to the external potential 
shows that the temperature gradient will cause a flow of particles; if they 
are charged, this will be interpreted as the thermo-electric effect. 

The rest of the paper is organised as follows. In §2 we outline the model 
in discrete space and time. It is a version of the Boltzmann equation, with 
discrete energy rather than discrete velocity. It is thus closer to |16] than 



[]17| . The collision operator is a bistochastic matrix T conserving energy and 
particle number (but no other quantities); it causes transitions (hopping) be- 
tween particles and holes which are nearest neighbours. The discrete system 
is thermodynamically consistent in its own right, and thus forms the natural 
discretisation of the continuum equations of motion eq. @. The state of the 
system at any time is described by giving the means of the 'slow variables', 
here taken to be the particle number, n = or n = 1, and the kinetic en- 
ergy, at each site. This information defines a unique grand canonical state 
at each site; the assumption that this is the state of the system after a small 
time-step is called LTE, the hypothesis of local thermodynamic equilibrium. 



It means that the state is specified by giving the density and temperature 
fields, which define a point on our information manifold, M. The dynamics 
of the state in one time-step is given by applying the map T, followed by 
projection of the resulting state back to ^A. 

In §3 the continuum limit of the dynamics is taken, with the help of 
MAPLE. The size of the lattice spacing i and the time step dt are arranged 
to satisfy dt = i"^, known as the diffusion limit. This ensures that the limit 
exists; it gives the dynamics above. 

In §4, we verify that the system can be put in Onsager form, though it 
is neither linear nor near equilibrium. Onsager symmetry is then seen to 
relate the Soret effect to the anomalous convection, which is thus the true 
dual or 'inverse' to the Soret effect. 

In §5 we summarise the paper, and conclude that our simple model, 
without any momentum or angular momentum, and with unrealistic density 
of states, is able to exhibit the Soret effect and its dual without the pain of 
solving a model with interparticle interactions. There is no reason to expect 
that these qualitative conclusions would be altered by a more elaborate 
model. For example, we could add a direct transfer of kinetic energy between 
occupied sites; this adds further terms to the energy current, and increases 
the diagonal part of the Onsager matrix, without affecting the cross terms 
(the Soret and Dufour terms). 

2 The Discrete Model with Hard Core 

We start with a finite lattice A C £Z^, where v is the dimension of space. 
A typical point of A will be denoted by x. At each site x there can be 
at most one molecule; this expresses the physical idea that the molecule 
has a repulsive core of diameter L We take a very simple model of the 
molecule; it has no spin or velocity, and its possible states are labelled by 
k G {0, 1, 2, . . .}, denoting the kinetic energy K, = fee; here e > represents a 
quantum of energy. This absence of velocity as a label for the point in phase 
space leads to a useful simplification of the mathematics compared with the 



discrete velocity Boltzmann equation [17|; we call it the Smoluchowski point 
of view. The sample space at x is thus 

O, = {0,0,1,2,...} = {0} UN. (6) 

Here, w^,. = £ ^x means that the site x is empty, and uox = fe^ G N means 
that the site x is occupied, and that its kinetic energy is kxt. We can argue 



that the fluctuation of the velocity of a real gas from its local mean value 
(the drift) is a fast variable, and is thermalised in one time-step. In this 
model, momentum is not accounted for; however, kinetic energy is. 

The sample space, also known as the phase space, for the system is taken 
to be 

n=l[n,; (7) 

the motivation for choosing this is given in [22|. The product structure 



of O ensures that the particles are indistinguishable, since a configuration 
is given by a field {uJx}x<=a- Thus the Gibbs paradox is avoided from the 
start, and there is no need to introduce symmetrising factors. In a classical 
stochastic description of such a system, an observable is a random variable, 
that is, a real function on Q. The set of all bounded random variables form 
an algebra, denoted here by A. At time i G N the state of the system 
is described by a probability measure on 0, which we denote by p; in the 
theory of C*— algebras, p is called a normal state on A. The set of normal 
states is a convex set denoted by S. 

A useful observable is the number of particles at x: 

. , . / if (Ja; = 0, , . 

^^(^^ = \ lif^^eN. ^^) 

The total number of particles is then the random variable 

A/-H = 5]Ar,H. (9) 

weA 

We introduce an external potential V{x) for the particles. Contrary to 
[^, we shall not find it necessary to introduce an interaction between the 
particles in order to get the Soret and Dufour effects. The full treatment of 
a model with interaction leads to a much more complicated theory than the 
present one. The 'mean field' treatment of the interaction can be handled 
as in Q; in this approximation, the interaction does contribute to the Soret 
effect, but is not the whole story. 

The total energy is the random variable 

£ = }C + V (10) 

where 

V(c^) = Y. VixWxicv); IC{co) = J2 ^kxMxiuj), (11) 



e > being an energy unit, and 

w = {uj^}x(zA, w^ = or k^. (12) 

We are interested in states with finite mean energy 

p. 6 := ''^p{u!)£{uj) < oo, (13) 

and finite partition function 

Z = Y, e-'^^(^) < oo for all /3 > 0; (14) 

this expresses thermodynamic stability. 

The function £ divides fi into disjoint energy shells: 

nE = W^^-- £{io) = E}, ^ = \J ^E- 

E 

Each VIe must be a finite set, because of eq. (p!4[). Similarly we can write 

^E = \_\^E,n-, 

n 

where 

^E,n = ^E n ^n, il„ = {w G fi : M{lo) = n} . 

The next step in constructing the model is to give a bistochastic map T 
on S, representing one time-step of the linear part of the dynamics. This 
replaces the 'collision term' of the full Boltzmann equation. In place of the 
large number of collision invariants of the Boltzmann equation, we require 
just two conserved quantities, the energy and the particle number. Now, T 
is determined by linearity from its action on the point measures 5^^^^', and 



these can be identified with the points of $7 |20|. To conserve £ and J\f, T 
must map each each il.E,n to itself. When we give T, we are specifying the 
conditional probability of transition from to to co'. We do not attempt to 
find the most general bistochastic map, but limit ourselves to the case where 

1- T^^,uj' = T^^',ui, the symmetric case. 

2. 7L,Lj' = if w and uj' differ at points x, x' € A which are not nearest 
neighbours. 



3. A transition occurs only by the movement of a particle; there is no 
direct transfer of kinetic energy between adjacent particles. 

To conserve energy, any change in potential energy in a transition must be 
balanced by the opposite change in kinetic energy. Since the latter is an 
integral multiple of e, we must suppose that all differences V{x) — V{y), 
with X, y G A, are integral multiples of e. We shall work out the model in 
detail when z^ = 1, and A is an 'interval' of contiguous points in iZ. Thus, 
for X G A with x + £ also in A, we can write 

V{x + i)- V{x) = ewx, with w^ G Z. (15) 

Because of (2) and (3), above, the possible transitions involve the movement 
of a particle from x to a hole at x + i oi x — i, or vice versa. We shall choose 
X, and consider the case Wx > 0; other cases are treated similarly. 

Suppose that there is a particle at x and a hole at x + i. In order for 
the transition x i— > x + £ to be possible, the kinetic energy of the particle at 
X must be at least ewx] after the transition its kinetic energy is €{kx — Wx), 
as the particle loses ewx in moving uphill. We call e{kx — Wx) the spare 
kinetic energy of the transition. If there is a particle at x + £ and a hole 
at X, the particle can move down to x whatever its kinetic energy happens 
to be. This move is the inverse to the first, and to arrive at a symmetric 
Markov matrix, we choose the transition rates for these two processes to be 
the same. We have argued |^^ that a good model for the transition matrix is 
to choose the rate to be proportional to the spare kinetic energy, by analogy 
with Einstein's hypothesis of stimulated emission. To this, we add one more 
unit (the 'spontaneous emission') which in the event makes no difference to 
the continuum limit. 

We see from our answer, eq. (^), that this choice leads to a linear increase 
in the thermal conductivity with temperature, roughly in accord with exper- 



iments [14|, p. 574, for various gases in the range 100° — 300° K. It also leads 
to a diffusion rate that increases linearly, ||2^, as chosen by Smoluchowski. 
Let X G A; let to and lo' be two sample points on the same f^_E,n, and let 
oj' = LOy for y / X or X + ^. Suppose that iOx = kx > Wx, and cux 



y 



+i 



ujx = and oj'x+i = kx — Wx- We call this the hopping criterion, and define 



the transition matrix Tx by 







00 U)' 


Tx = 


-- UJ 


1 - 2\e{kx -Wx + l) 2Xe{kx - 




J 


2Xe{kx -Wx + l) 1 - 2Xe{kx - 



Wx + 1) (16) 

Wx + 1) 



It was noticed in [E^, 24 1 that transition rates which grow with energy need 



care when the time is discrete, as the transition probabihty becomes larger 
than 1 for large energies; the stay-as-you-were probability can become neg- 
ative. Not only is this unphysical; it leads to unstable difference equations 
and spurious chaos ||l^. It is therefore desirable in our model to cut off the 
hopping probability to zero above 

km.. = (2Ae)-i - 1. (17) 

By making e smaller we can make this cut-off as large as we please. For 
simplicity of notation, we put the matrix elements of T^ equal to zero unless 
CO and to' obey the hopping criterion. That is, we remove the diagonal line 
of units. A matrix like T^ can be defined for each x G A. In particular, T^-i 
describes the move down from x to x — i, and the move up from x — i to x, 
which is assigned the probability 

2Xe{k^_i - w^_e + 1) = 2Xe{k^ + 1), (18) 

or zero if k^-e < Wx-i- 

Both Tx and Tx-n affect the population at x, but other transitions, 
Ty, y ^ x,x — i, do not. We define the linear part of the local dynam- 
ics at x to be {Tx + Tx-e)/2, and the full linear dynamics to be given by the 
symmetric Markov matrix 

T = Y.{Tx + Tx-i)/2. (19) 

The next step in statistical dynamics, the analogue of the Stosszahlansatz, 
is the specification of the thermalising map Q. This projects the result 
of one time-step, p i— > Tp, onto a subset of states called the information 
manifold. We seek to implement mathematically the physical idea that 
the kinetic energy at a point x completely thermalises in one time step. 
The description of the local state by the full distribution of the random 
variable K,x is replaced by one variable, the temperature 0(x). According 
to ll^, |ll we have to choose a subspace X of slow variables, spanned by 
{Xq = 1, Xi = H, X2, . . . , Xj\f}, in the notation of |2^. Given a state p, we 
record the mean values 

Vj=P-Xj, j = l,2,...,N; (20) 

we then define Qp to be the state of maximum entropy having these mean 
values. The set of states A4(Af) = {Qp : p € S} has a Riemannian structure. 



and is called the information manifold of A'. It is parametrised by the means, 
which are called the mixture coordinates. If we start with p G M, we define 
one time-step of the full (non-linear) dynamics to be p i-^ QTp. This maps 
A4 to itself, and is given by difference equations in N variables. This is far 
fewer than needed for the linear dynamics. 

In the present model, we take the slow variables to be the span of {A^ : 
X E A} U {ICx : a; E A}, a total of A^ = 2|A| variables. In particular the total 
energy is a slow variable. Let us denote the means of Afx and JCx by the 
fields Ux and Kx'- 

nx=p.Afx] Kx=pXx. (21) 

The algebra of observables. A, is the tensor product 

A = ^xAx (22) 

of local algebras. The state Qp is then the state of maximum entropy with 
these means; it is well known that this is the state in LTE, the product of 
local grand canonical states, independent over A; thus 

Qp{L0) = l[px{i0x) (23) 

X 

where Px is the marginal probability of p on i^x- In its turn, px is the state 
on Qx of maximum entropy among states with given means Ux, Kx- It is 
therefore the (fermionic) grand canonical state, and so has the form, for 
some parameter Px'- 

Pxm = l-nx (24) 

Px{k) = UxZ-^e-'^^-, (25) 

where Zx = {)- — exp(— e/?a;))~^. We determine (3x or the temperature Qx = 
l/^ksPx), from the mean-value Kx'- 

K{x) = p-lCx = enxY.kZ-^e~'^^- (26) 

fe=0 

eux 



o€f3x 



(27) 



Here we have taken the sums up to infinity, neglecting the term exp(— /JAjmaxc) 
compared to 1. 



We note that 

n{x) = ^p.(fe) (28) 



fc>0 



e 



-€/3 



{l + en/K)-^ (29) 

Z = l + K/{en). (30) 



Our strategy for specifying one time-step in the dynamics is to start 
with p ^ M., and so of the form eq. (|2^), (p^), (^0|), and then find p = Tp, 
or rather, its marginals. Prom the marginals of p we find the new values 
ri(x), K{x) of the mixture coordinates; from these we can find the new 
exp(e/3) = 1 + en/K. It remains to explain how to find the marginals of j5. 

For a product state over A, eq. (|23|) , the probability of the set of points 
like w, satisfying the hopping criteria at x, (i.e. x occupied with kinetic 
energy ek^, and x + l empty) is p{x, k){l — rix+i). Similarly, the set of points 
like a;', (i.e. with x empty, and x+i occupied with kinetic energy e{kx—Wx)) 
is p{x + i,kx — tfx)(l — nx)- The other sites y ^ x, x + i are unaffected by 
Tx and these factors in eq. ( p3|) can be summed over ujy to give unity. Thus 
the marginal at x of TxP is 

Txp{x,k) = p{x,k) -p{x,k){l-nx^£)2Xe{k + l-Wx) 

+ p{x + i,k-Wx)il-nx)2Xe{k + l-Wx) ■ (31) 

Similarly 

Tx-ep{x,k) = p{x,k) -p{x,k){l-nx-e)2Xe{kx-e-'Wx-e + l) 

+ p{x - i, kx-e) (1 - Ux) 2A {kx-e - Wx-e + 1) . (32) 

Here, k = kx-i — Wx-t must hold, to conserve energy in the transition. So 
the change in the marginal at x due to {Tx + Tx-i) /2, and therefore to T, 
simplifies a bit to 

p{x,k) = p{x,k) -p{x,k){l-nx-e)^e{k + l -Wx) 

+ p{x + i,k — Wx) (1 — Hx) Ae (A; + 1 — Wx) 

- p{x,k){l-nx-e)Xe{k + l) 

+ p{x -e,k + Wx-i){l-nx)Xe{k + l). (33) 

From this we can find the new values, the one-step updates, n, K: 

nx = ^pix,k) 

k>0 



10 



n.x - Aen^. I i - n^+H z^' > [K + i-w^)e '''^" 



ix - Aen^ (1 - n^+i) Z^^ ^ (fe + 1 



k>Wx 



Wx e 



-e{k-Wa:)l3x + i 



+ Aen^+<? (1 - rix) Z^^^ ^ (/c + 1 

k>Wa: 

- Xeux (1 - n,_,) Z-i ^ (A; + 1)6"^^^- 

+ Aen,_,(l-n,)Z-i,5](fc+l)e-^(^+--^)'^-^ (34) 

fe>0 



and 



Kx = Yl (^kp{x, k) 

k>0 

= Kx- e^Xux (1 - n^+i) Z'^ ^ k{k + l-Wx) e-'''^^ 
+ e^Xrix+i (1 - n,) Z'^/, ^ fc (fc + 1 - «;,) e-^^^— -)^/5-+^ 
- e^An, (1 - rix-e) Z'^ ^ k{k + l)e~'^^^ 

k>0 

+ e^An,., (1 - n,) Z-i, ^ fc(fc + l)e-^('=+"'-^)^-^ (35) 

A:>0 

The sum should go up to the largest k consistent with the positivity of the 
diagonal matrix element of T^j^^^,/ at x and x — i. This ensures that the model 
is entropy increasing, and that the discrete dynamics is stable. 

The dynamics is the explicit map n >-^ n, K >-^ K, and from K we can 
compute P, and thus get the new point of A4. In the next section we take 
the continuum limit of this dynamics. The numerical solution of the re- 
sulting reaction-diffusion equations might be best done using these discrete 
equations: energy is conserved and entropy increases in all approximations, 
which ensures that the numerical solutions have reasonable physical proper- 
ties, and do not exhibit spurious chaos, (chaos not present in the continuum 
equations). 

3 The Continuum Limit 

We now take the continuum limit of the updating equations, thus: 

dt = f ^0; e = -ii^O (7 fixed) (36) 

11 



K -K 

Jt 

V{x + (.)-V{x) 



dK 

dt ' 


1 

J -^ Pm', 


n 


dV 
dx 


and so on. 





(37) 
(38) 



This is achieved by writing 

p,±i = P ± 1(31 + (1/2) f 132 (39) 

n^±e = n±inl + {l/2)e^n2 (40) 

ew^-e = V{x) - V{x - £) = evi - {l/2)i^V2 (41) 

h-n = A1 + A2 + A3 + A4. (42) 

Here, pi is first the derivative of (3, 02 is the second derivative, and the 
same for the functions n and V; the A's are four expressions in eq. (p^): 

Al = -eAn(l-n3)e-^'"^(l-e-^^)^(A: + l)e-^^'= (43) 

fc>0 

A2 = eAn3(l-n)(l-e-^^=^)^(A; + l)e-^^=''= (44) 

fc>0 

^3 = -eAn(l - n4) (l - e-'f^) ^ ('t + 1)6""^*^ (45) 

fc>0 

AA = eAn4(l-n)(l-e-^^^)e-^^^"'^^(A: + l)e-^^^^ (46) 

fe>0 

where /33 = Px+i^ /?4 = Px-e, n3 = Ux+e, nA = Ux-e, wA = Wx-e- The sum 
is ahowed to go to infinity in the continuum limit. 

We ask MAPLE to evaluate n — n to lowest non-vanishing order in i, 
namely, £'^; we substitute e = 7^, with 7 fixed, n3 = n + in' + (l/2)^^n" 
etc., P(x) = 1/0(2;), and n — n by d^dn/dt. In one dimension the result, as 
^ ^ 0, is 

OTi 

— + div> = (47) 

where 

jn = -A (Gn' + (G' + V) n(l - n)) . (48) 

We see that this is independent of 7. 

The particle density is n/i. = p; then n = pi which we replace by p/pm 
rather than by zero in the limit. So by dividing eq. (^) by i gives our 
equation of motion for the particle density: 

jc = -A (GVp + p{l - p/pm)V{e + V)) . (49) 

12 



The term Xp{l—p/pm)V@ is the thermal diffusion; if p <C Pm we see that the 
Soret coefficient is exactly A, but that for dense fluids (near solidiflcation) it 
is smaller. The Soret term is absent in the model of dense fluids presented in 
[p2| . This arises because there the particles do not carry heat, unlike those 
of the present model. 

For the kinetic energy density K{x), we write 

k-K = AI + A2 + A3 + A4:, (50) 

where 

Al = -e'^Xn{l-n3){l-e-"^)e-''^'"Y,{k + w){k + l)e-'^'' (51) 

fe>0 

A2 = e2A(l-n)n3(l-e-^^3^^(^^^)(^^l)g-e/33fc ^52) 

fc>0 

A3 = -e2An(l - n4) (1 - e"^^) ^ A:(A; + l)e-^^^ (53) 

fc>0 

Ai = e2An4(l-n)(l-e-^^^)e-^^^"'^^A:(A; + l)e-^^". (54) 

fc>0 

Now we put K = n@ and take the limit £ — > 0; MAPLE gives 



d{ne) 
dt 



-A <^ n'^iV'f - SV'n'e - 3V'ne' - n{V'f 

+ An'^ie'f + An^ee" - 4nee" + 2t?QV" 
+ 3nVG' + Snee'n' + 4nn'y'e 

- 2n"G2 - Sn'ee' - 2nQV" - 4n (G') 

= -jnV - 2div (ei„) + 2Adiv (6(1 - n)n<d)') . (55) 

Again we see that 7 drops out. Putting n/i = p and l/£ = pm we get in 
i/-dimensions: 

^^ = -Jc'^V - 2div (Gje) + 2Adiv (Gp(l - p/p„,)Ve) . (56) 

ot 

Thus, at least at a formal level, the equations discussed in §1 are 'derived' 
from a discrete stochastic model. Our derivation has not proved that the 
solutions to the discrete equations converge to solutions of the nonlinear 
coupled equations in the continuum. Indeed these are not uniformly elliptic, 
and a proof of existence of even local solutions needs care. Some results for 



similar systems are presented in |25 
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4 The Onsager Form 

We now argue that only the 'anomalous' part of the Dufour effect is the true 
Onsager dual to the Soret effect. We see this in the model of [22|; this model 
has no Soret effect, but the heat current does contain the normal convection 
term @jc- In the classical Boltzmann equation the heat content of a fluid is 
3/20/>, instead of our Qp; we would expect, then, a convection term 3/20 jc 
in the heat current. However, in Q the term 5/20jc is found. We see that 
of this, 3/20jc is normal convection, and Qjc is anomalous convection; the 
latter is the same as in the present model. We shall show that Onsager 
symmetry relates the Soret effect to the anomalous part of the convection. 
Both the model of [^], and the present model, have the same sample 
space, slow variables and conserved quantities, namely the particle number 
and energy. It follows from the argument given in |^3| that both models have 
the same entropy and thermodynamic forces. They differ in their dynamics, 
which in Onsager theory is determined by the way the currents depend on 
the forces. We first find the entropy of the discrete model, and then take its 
formal continuum limit. The entropy of the state 

P(^) = ni(l - -^-(^))(l - n,) + M,{u;)p,{u;,)} (57) 

X 

which is independent over A, is the sum of the contributions at each x; thus 

UJ 

= -X!(^ -na;)log(l -n^) - ^Hx^ Px{k) log pxik) 
= -^{l-nx)log{l-nx) -^rixlogrix 

X X 

+ ^n,,{logZ;3.+/3xe/(e^^^-l)} (58) 

X 

Now, 

-^rixlogrix - ^(1 -na,)log(l -n^.) (59) 

X X 

and in the continuum limit, ^ — > 0, X^x ^ becomes / dx, Uxji becomes p{x) 
and \ji becomes pm- So the dichotomic part of the entropy becomes the 
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differential entropy 

p{x) log p{x) dx - I {pm - p) log (1 - p/Pm) dx, (61) 

apart from the large positive term 

-logi I p{x)dx. (62) 



However, this term can be dropped; for / p{x) dx is constant in time, so the 
divergent term does not contribute to dS/dt. The Gibbsian term in Sx also 
has a simple limit: 

nx{\ogZp^+(3xKx) = n,(-log(l-e-^^-)+/?,e/(e^'^-l)) 

= -n^^loge-n^. log /3j, + na; + 0(e). (63) 

Again, summing over x leads to the divergent but constant terms (— log e + 
1) X^x "^x: which can be dropped, leaving 

- ^ log PxUx -^ / p{x) log 0(x) dx; 
The continuum entropy is therefore 

S = - p{x) log p{x) dx- j {pm - P) log (1 - p/Pm) dx 

p\ogQ{x) dx, 



as claimed in |g2[, eq. (49). 

The next step |23] in setting up a comparison with Onsager theory is to 
write 5 as a (continuous version) of Onsager's ansatz [^] 

e5 = ^X"j„. (64) 

a 

We use the identities 

j pQ/Qdx = [{pe)/Gdx (65) 

pp/pdx = (pm - P) P/ (PmiPm - P)) dx = (66) 
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to get 

S = - plogpdx+ plog{l- p/pm) dx 

+ I p\ogQdx+ I {pQ)/Qdx. (67) 



In both the present model and that of [22| there are two conserved densities, 



the particle number p and the energy £ = p{V + 0); thus pG = £ — pV . 
We impose the condition of 'no flow' on the boundary, 9A, which could be 
at infinity. Thus the components of the current normal to the boundary, j^ 
and jjf, vanish on dK. In both models, then, we use the conservation laws 
and then integrate by parts in eq. (^) and discard the boundary term, to 
get: 

S = - j p{\og{ppJ{p^-p))-\oge) + j (£-pV)/@dx (68) 

= - j JcV (log(p/(l - p/pm)) - loge + V/Q) dx 

+ I ie.V{l/Q)dx. (69) 



This is true whatever the equation of motion, provided that the particle 
number and energy are conserved. In particular, we can vary the hopping 
probability. A, as a function of (x, t), so that jc and jg are arbitrary functions. 
So from the Onsager ansatz for S we can read off the thermodynamic forces 

1 



p{l- p/pm) 



These expressions for the thermodynamic forces conjugate to the currents 
jc and je hold for both models; we consider them in turn. In the model of 
[P2| the equations of motion are 

Jc = -x{vp + pil-p/pm)^^ (73) 

je = -X{p{l-p/pm)^e} + {e + V)jc. (74) 
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This model has no Soret effect, but shows the 'normal' convection, in our 
terminology, because of the contribution Qjc to the heat current. This 
will be observed as the Dufour effect, since there is a heat flow if VG = 
but V/9 ^ 0. Somewhat luckily (for Onsager theory), the currents are 
linear expressions in X^ and X^, with no derivatives, but with nonlinear 
coefficients: 

je = Xp{l-p/pm)^ + XiV + e)pil-p/pm)^; (75) 

je = {V + e)Xp{l-p/pm)^ 

+ X{Xp{l-p/pm)@^ + {e + Vfp{l-p/p^m)}^. (76) 

We see that Onsager symmetry holds, and that the Onsager matrix is posi- 
tive definite. Hence entropy is an increasing function of time, as expected. 
In the present model, we see that the currents, given in eq. (|2|), can be 
expressed in terms of the same thermodynamic forces eq. (|72|) thus: 

jc = Xp{l-p/pm)^ + \{V + 2@)p{l-p/p^)^ (77) 

j, = ^v + 2e)Xpil-p/p^)^ 

+ xUv + 2efp{l-p/pm) + 2e^p{l-p/pm)^\. (78) 

Again, Onsager symmetry and positivity hold. The anomalous part of the 
convection, the factor '2' in y + 20, can be inferred using Onsager symmetry 
and the Soret term, VG, in jc- In our approach, all these properties, the 
Soret effect, the anomalous Dufour effect and the Onsager symmetry, follow 
from the model, rather than being put in, as in the Onsager theory. 

We can relate the thermodynamic forces to the gradients of the canonical 
coordinates ^"(x) of the information manifold (when the states are indepen- 



dent over A), according to the general theory; see eq. (15) of p3|]. We shall 
now verify that this is true in this model, in the continuum limit. Recall 
that in the discrete model we write the density matrix as 



P 



Eo^"^'^. (79) 



Here, Hq = 1, ^'^ = logH, where H is the grand partition function; the Ha 
are the conserved densities. In our model, these are the number of particles 
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at X, Mxi and the energy aX x, £x = Vx + ^x- A point on the information 
manifold, for finite A, is the product over x of a probabihty which we write 
in the usual form of a local equilibrium state: 

p(x) = e'^°s^e-'^-(^--^-^-\ (80) 

We can thus identify the canonical coordinates as ^% = —(3xfix and, as 
expected, ^| = (3x- To relate ^^ to n, one of the mixture coordinates, note 
that (omitting the label x) 

E-^ = l-n (81) 

and 

E = Y.ef'^^e-"' (82) 

u> 

k>0 

e/3(M-v) 

= i + r^i^- (83) 



We combine eq. ( pl|) with eq. (|83|) to get 

p^ = /3y + log (n/(l - n)) + log (l - e"'^) (84) 



= pV + \og ^, — + loge-log/? + 0(e) (85) 

P - P/Pm 

for small e. In the continuum limit we drop the infinite constant loge, as 
only the gradient is used. Thus we get 

r = -/?// = -log ^ +loge-^ (86) 

1 - P/Pm B 

and X'^/Q = V^'^, using eq. (70). This, together with eq. (72) shows that 
(CSC") are "potentials" for {-Xye,-X''/e). 



5 Conclusion 

We have constructed an example of nonequilibrium thermodynamics obey- 
ing the first and second laws, and which exhibits the Soret and Dufour 
effects. Apart from the hard core, no interparticle potential is postulated, 
and indeed the effects persist in the limit pm -^ cxd, corresponding to no 
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hard core. This should be compared with the kinetic theory described in 
[^, p. 103; this gives the impression that the Soret and Dufour effects depend 
on the careful inclusion of interparticle forces and that they are present only 
in gas mixtures. In our approach the effect of the interparticle forces are 
included only indirectly, in that the dynamics includes the thermalising map 
Q, which ensures that the motion is confined to the information manifold. 
Thus, whatever the forces are, they keep the system in local thermodynamic 
equilibrium. This leads to a simple understanding of the effects; regions of 
higher temperature contain more high-speed molecules than regions of low 
temperature, and so they preferentially move from high to low tempera- 
tures. In our model there is no velocity variable, and by high speed we 
mean molecules that hop with greater probability. In Q an abnormal con- 
vection is also found, and is interpreted as "convection of enthalpy" . It is 
the differential Soret effect between isotopes that is emphasised there. This 
is the reason for the great industrial importance of the effect; it allows gas 
mixtures to be separated by a heat gradient. To get this effect we would 
need to allow the hopping rate A to depend on the mass of the molecule. 

We have remarked that the density of states in our model is unity; in a 
semiclassical model in one dimension, there are two states of each energy, 
corresponding to the two directions of motion, so apart from this trivial 
factor, it would appear that our model is one dimensional. However, the 
hopping rate is not that of a particle in one dimension. In the free dynamics 
between collisions the rate of movement is proportional to the velocity i.e. 
E^'"^ , whereas in our model it is proportional to the spare kinetic energy E. 
This larger hopping rate partially compensates for the lack of multiplicity 
of states, as can be seen by the following argument. In three dimensions the 
number of states of a single particle in the semiclassical model is propor- 
tional to the volume of phase space, and thus has a factor p^ dp (where p is 
now the momentum) which is proportional to E^''^dE. Suppose a particle 
in three dimensions at a lattice site x G A having momentum p hops to a 
neighbouring site, x + i against the force given by the gradient of the poten- 
tial V. Its momentum after the transition will be slightly different, say p', 
which is determined by p and the force. Thus each of the many states at x 
can make a momentum-conserving transition to only one of the many states 
of the right energy at x + i. So the number of states making the transition is 
proportional to the number of states, namely E^''^. The rate of flow of such 
an element of phase space is proportional to E^''^ as well, so the number 
making the transition is En{E), where n{E) is the occupation number. This 
is the same as that due to stimulated hopping, as assumed in the present 
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paper in eq. (18). As we saw, the extra hop due to spontaneous hopping 



did not contribute to the continuum hmit. It should be mentioned that 
this does not mean that model of the present paper gives the same answer 
(apart from trivial factors) as a three-dimensional model with the "correct" 
multiplicity, and rates proportional to the speed. Indeed, the factors Z~^ in 



the rate equations eq. (|3J) and (|35D also depend on the multiplicity. These 
factors are proportional to 0~'^/^ instead of 0~^ as in the present paper. 

We can also regard the present model as describing a dense liquid of com- 
plex atoms, which have little kinetic energy, but which have a large number 
of excited states modelled by the levels e/c. A; = 0, 1, 2, ... , fcmax- As re- 
marked, to get the equations (1) and (2), we approximate (1— exp(— /Je/cmax)) 
by 1, which requires that the temperature is low enough so that the states 
of high energy are not excited much. The density of states can be indepen- 
dently checked by experiment. For such a liquid, it would be interesting to 
see whether it is true that Dufour effect is double the convection (the factor 
2 in eq. (2)). 

Our model is more in the spirit of the discrete energy Boltzmann equa- 



tions [^] than the discrete velocity models [17| which conserve momentum 
as well as energy. For example, pi considers a lattice gas version of the 
Boltzmann equation, with a collision term between pairs of particles. The 
Dufour effect is predicted. In |p, only two speeds occur. The number of 
particles having each speed is separately conserved, so there is no thermal 
mixing in the scattering. The up-date equations for the densities are of 
the eighth degree. The authors compare numerical simulations of the exact 
model with solutions to a simplified model in which the momentum is put 
equal to zero, and a thermalising assumption is imposed. They remark that 
the simplified model is very close to the simulations. Our result shows that 
the Soret effect as well as the Dufour effect can be obtained without the 
two-body scattering. 

We may rather easily vary the lattice shape, and allow hopping to next- 
nearest neighbours; we tried several such variants, and got the same contin- 
uum limit provided that we adjusted A so that the finite-difference operator 
in the transition matrix T approximates the Laplacian. So the limit is rather 
robust. 
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6 The Maple V Program 

6.1 The Density 

tJWe use the notation nl = dn/dx, n2 = d^n/dx^, 

tj n3=n(x+l), n4=n(x-l), and the same for 

tj beta and w. 

assume(epsilon>0);assume(beta>0); 

Al : =-lambda*epsilon*n* ( l-n3) *exp(-epsilon*beta*w) * (1-exp (-epsilon*beta) ) * 

sum((k+l)*exp(-epsilon*beta*k),k=0.. infinity); 

Bl:=simpnfy(Al); 

Cl:=series(Bl,epsilon,5); 

Dl:=convert(Cl,polynoni); 

Fl:=simphfy(Dl); 

Gl:=subs(epsilon=l*gamma,Fl); 

Hl:=subs(n3=n+l*nl+(l/2)*r 2*n2,Gl); 

Jl:=series(Hl,l,3); 

Z 1 : =convert ( J 1 ,polynoni) ; 

assume (b et a3 > ) ; 

A2: =epsilon*lambda*n3* ( 1-n) * ( l-exp(-beta3*epsifon) ) *sum( (k+ 1) * 

exp(-beta3*epsifon*k),k=0.. infinity); 

B2:=series(A2,epsilon,5); 

C2:=convert(B2,polynom); 

D2:=simphfy(C2); 

F2:=subs(epsilon=l*gamma,D2); 

G2:=subs(n3=n+l*nl+(l/2)*1^2 *n2,F2); 

H2:=subs(beta3=beta+l*betal+(l/2)*r2*beta2,G2); 

J2:=series(H2,l,3); 

K2:=convert(J2,polynom); 

Z2:=simplify(K2); 

A3: =-epsifon*lambda*n* ( l-n4) * ( 1-exp (-beta*epsilon) ) *sum ( (k+ 1 ) * 

exp (-beta*epsifon*k) ,k=0. .infinity) ; 

B3:=series(A3,epsilon,5); 

C3:=convert(B3,polynom); 

D3:=subs(epsilon=l*gamma,C3); 

F3:=subs(n4=n-l*nl+(l/2)*r *n2,D3); 

G3:=series(F3,l,3); 

H3:=convert(G3,polynom); 

Z3:=simplify(H3); 



23 



assume (b et a4 > ) ; 

A4: =epsilon*lambda*n4* ( 1-n) * ( l-exp(-beta4*epsilon) ) * 

exp(-beta4*epsilon*w4)*sum((k+l)*exp(-beta4*epsiloii*k),k=0.. infinity); 

B4:=series(A4,epsilon,5); 

C4:=convert(B4,polynoni); 

D4:=subs(epsilon=gamma*l,C4); 

F4:=subs(n4=beta-l*nl+(l/2)r 2*n2,D4); 

G4:=subs(beta4=beta-l*betal+(l/2)*r2*beta2,F4); 

H4:=subs(w4=w-l*wl+(l/2)*r2*w2,G4); 

J4:=series(H4,l,3); 

K4:=convert(J4,polynoni); 

Z4:=simplify(K4); 

ansl:=simplify(Zl+Z2+Z3+Z4); 

ans2:=subs(beta=l/Theta(x),ansl); 

ans3:=subs(betal=diff(l/Theta(x),x),ans2); 

ans4:=subs(beta2=diff(l/Theta(x),x$2),ans3); 

ans5: =subs(w=diff (V(x) ,x) ,ans4) ; 

ans6:=subs(wl=diff(V(x),x$2),ans5); 

ans:=siniplify(ans6); 

Simple manipulation then gives eq. (|l]). 

6.2 The Heat 

tj Soret2.ms; as before, we take the four terms from A to Z. 

tJ betal denotes dfi/dx, beta2 denotes (f'P/dx^ , 

jl beta3 denotes j3{x + I)] beta4 denotes (3{x — i). 

H and the same for n and w. 

assume(beta>0) ; assume(epsilon>0) ; 

Al:=-lambda*epsilon^ 2*n*(l-n3)*exp(-beta*epsilon *w)* 

(l-exp(-beta*epsilon))*sum((k+w)*(k+l)*exp(-beta*epsilon*k),k=0.. infinity); 

Bl:=simplify(Al); 

Cl:=series(Bl,epsilon,6); 

Dl:=convert(Cl,polynom); 

Fl:=subs(epsilon=l*gamma,Dl); 

Gl:=subs(n3=n+l*nl+(l/2)*r2*n2,Fl); 

Hl:=series(Gl,l,3); 

Jl:=convert(Hl,polynom); 

Zl:=simplify(Jl); 

assume (b et a3 > ) ; 
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A2:=epsilon^2*beta*(l-n)*n3*(l-exp(-epsilon*beta3))* sum((k+w)*(k+l)*exp( 

epsilon*beta3*k) ,k=0. .infinity) ; 

B2:=simplify(A2); 

C2:=series(B2,epsilon,6); 

D2:=convert(C2,polynom); 

F2:=subs(epsilon=l*gamma,D2); 

G2:=subs(n3=n+l*nl+(l/2)*r2*n2,F2); 

H2:=subs(beta3=beta+l*betal+(l/2)*r2*beta2,G2); 

J2:=series(H2,l,3); 

K2:=convert(J2,polynoni); 

Z2:=simplify(K2); 

A3:=-epsilon^2*lambda*n*(l-n4)*(l-exp(-epsilon*beta))*sum(k*(k+l)* 

exp (-epsilon*beta*k) ,k=0. .infinity) ; 

B3:=simplify(A3); 

C3:=series(B3,epsilon,6); 

D3:=convert(C3,polynom); 

F3:=subs(epsilon=l*gamma,D3); 

G3:=subs(n4=n-l*nl+(l/2)*1^2*n2); 

H3:=series(G3,l,3); 

K3:=convert(H3,polynoni); 

Z3:=simplify(K3); 

assume (b et a4 > ) ; 

A4: =epsilon ^ 2*lambda*n4* ( 1-n) * ( 1-exp (-epsilon*beta4) ) * 

exp(-epsilon*beta4*w4)*sum(k*(k+l)*exp(-epsilon*beta4*k),k=0.. infinity); 

B4:=simplify(A4); 

C4:=series(B4,epsilon,6); 

D4:=convert(C4,polynoni); 

F4:=subs(n4=n-l*n2+(l/2)*1^2*n2,D4); 

G4:=subs(beta4=beta-l*betal+(l/2)*1^2*beta2,F4); 

H4: =subs (epsilon=l*ganinia, G4) ; 

J4:=subs(w4=w-l*wl+(l/2)*1^2*w2,H4); 

K4:=series(J4,l,3); 

L4:=convert(L4,polynom); 

Z4:=simplify(L4); 

ansl:=simplify(Zl+Z2+Z3+Z4); 

ans2:=subs(w=Vl/gamma,ansl); 

ans3:=subs(wl=V2/ganinia,ans2); 

tj Note that tlie result is of second degree in 1 and is independent 

tj of gamma. 
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ans4:=subs(beta=l/Theta(x),ans3); 

ans5:=subs(betal=diff(l/Theta(x),x),ans4); 

ans6:=subs(beta2=diff(l/Theta(x),x$2),ans5); 

Ans:=simplify(ans6); 

This gives us eq. (p5[). 
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